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Abstract 

In this note we prove NP-hardness of the following problem: Given 
a set of matrices, is there a convex combination of those that is a 
nonsingular M-matrix? Via known characterizations of M-matrices, 
our result establishes NP-hardness of several fundamental problems 
in systems analysis and control, such as testing the instability of an 
uncertain dynamical system, and minimizing the spectral radius of an 
affine matrix function. 



1 Introduction 

Several problems in the biological, physical, social sciences and engineering, 
can be modelled via matrices that possess special structure. A commonly en- 
countered case involves the use of nonsingular M-matrices, which are positive 
stable matrices with all their off-diagonal entries nonpositive. M-matrices 
enjoy a remarkably rich theory, with connections to the Perron-Frobenius 
theory of nonnegative matrices, and with many applications in different 
scientific fields. The books of Horn and Johnson (1991) and Berman and 
Plemmons (1994) provide extensive treatments of the subject. 

Negated M-matrices constitute a special case of Hurwitz matrices that 
play a pivotal role in the stability of dynamical systems (Khalil, 2002). A 
fundamental problem in this area is the analysis of the stability of a system 
under model uncertainty, with important applications in systems engineering 
and biology (Doyle and Csete, 2011). Often this uncertainty manifests as a 
polytope of system matrices, whose overall stability or instability needs to 
be established (Chesi, 2010). In this context, it is of interest to know what 
is the computational 'effort' that it takes to test the (in)stability of a given 
uncertain system. 



1 



In this note we prove NP-hardness of the problem of testing whether 
there exists a convex combination of a given set of matrices that is a nonsin- 
gular M-matrix. Using known characterizations of M-matrices, our result es- 
tablishes NP-hardness of several other problems with practical significance, 
such as the problem of testing the instability of a linear time-invariant sys- 
tem with polytopic parameter uncertainty, and the problem of minimizing 
the spectral radius of an affine matrix function. NP-hardness of the latter 
problem was independently established recently (Fercoq, 2011). 

2 The POLYTOPE-M-MATRIX problem 

A real matrix whose off-diagonal entries are all nonpositive is called a Z- 
matrix. If a Z-matrix is positive stable, that is, all its eigenvalues have 
positive real part, then it is called a nonsingular M-matrix (Berman and 
Plemmons, 1994, chapter 6). 

In this section we prove NP-hardness of the following decision problem: 

Problem 1 (The polytope-M-matrix problem). Given a finite set of 
real matrices, is there a convex combination of those that is a nonsingular 
M-matrix? 

Theorem 1. The polytope-M-matrix problem is NP-hard. 

Proof. We establish a polynomial-time reduction from the independent- 
set problem. This problem asks, for a given undirected graph Q = (V, E) 
and a positive integer j < \V\, whether Q contains an independent set V 
(a set of pairwise non-adjacent vertices) with size \V'\ > j. This problem is 
NP-complete (Garey and Johnson, 1979). 

An instance of polytope-M-matrix takes as input k real nxn matrices 
Aj, with i = l,...,k, and asks whether there exists a nonnegative vector 
7T = (tti, . . . , 7Tfc) T with Yli=i n i = 1> such that the matrix = J2i=i ^i-^-i 
is a nonsingular M-matrix. We will prove NP-hardness of the problem for 
the special case in which k = n, that is, when the number of input matrices 
equals their dimension. For any input instance of INDEPENDENT-SET , our 
reduction will construct an instance of polytope-M-matrix , such that a 
polynomial-time algorithm for deciding the latter would imply a polynomial- 
time for deciding the former. 

Let (Q, j) be an instance of independent-set , and let C be the nxn 
adjacency matrix of the graph Q. The matrix C is a symmetric zero-one 
matrix with all zeros in the main diagonal. Let Cj denote the i'th column of 
C, and let e^ denote the length- n vector with 1 in the i'th entry and all other 
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entries zero. The reduction constructs n block matrices Aj, for i = 1, . . . , n, 
of size (n + 1) x (n + 1): 



A 



I -(ej + Cj 



4 — _ P T 1 



(1) 



where I is the n x n identity matrix. Note that each Aj is a Z-matrix. 

Suppose there exists a polynomial-time algorithm that allows us to find a 
vector 7r = (7Ti, . . . , 7r n ) T in the probability simplex A, for which the matrix 
B,,- = Y^=i ^iA-i is a nonsingular M-matrix. Using the definition of Aj from 
(1), the matrix B,,- can be written as 



I -(I + C)7T 
7T 



'7T ^.J 1 

3 



(2) 



A necessary and sufficient condition for a Z-matrix to be a nonsingular M- 
matrix is that all its leading principal minors are positive (Berman and 
Plemmons, 1994, p. 135, condition E17). Consequently, since the (1, 1) block 
of B„- is the identity matrix, the matrix B„- will be a nonsingular M-matrix 
if and only if det(B,r) > 0. Using the Schur determinant formula for block 
matrices (Berman and Plemmons, 1994, exercise 6.1), the determinant of 
B,,. can be expressed as 

det(B 7r ) = i-7rT(I + C)7r, (3) 

and hence det(B,r) > implies 7r T (I + C)tv < i. Moreover, for any graph Q 
with adjacency matrix C, the following is true (Motzkin and Straus, 1965): 

miny T (I + C)y, (4) 



a{Q) yeA 

where a{Q) is the size of the maximum independent set of Q. Using (3) 
and (4), we conclude that the existence of a vector tt 6 A that satisfies 
det(B 7r ) > implies < 4, or a(Q) > j, and hence it implies the 

existence of an independent set V C V in Q with size \V'\ > j. 

Since all steps in the above reduction take time at most polynomial in the 
size of the problem input, we conclude that the existence of a polynomial- 
time algorithm for deciding polytope-M-matrix would also allow us to 
decide independent-set in polynomial-time, implying P=NP. □ 
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Remarks: Effectively the proof hinges on the fact that nonsingular M- 
matrices are closed under Schur complementation (Horn and Johnson, 1991, 
p. 128, problem 15), a property that they share with the class of positive 
definite matrices. However, unlike positive definite matrices, M-matrices 
are not closed under regular matrix addition, accounting for the difficulty of 
optimization problems involving M-matrices. If the input matrices are sym- 
metric, the problem is convex, since symmetry of a nonsingular M-matrix 
implies its positive definiteness (Berman and Plemmons, 1994, p. 141, ex- 
ercise 2.6), in which case the corresponding optimization problem can be 
solved by semidefinite programming (Boyd and Vandenberghe, 2004). 

3 Corollaries 

Berman and Plemmons (1994) provide a list of more than 50 equivalent con- 
ditions for a Z-matrix to be a nonsingular M-matrix. Our result establishes 
NP-hardness of all 'polytopic' versions of them. For instance, using condi- 
tions Di6 and N38 of Berman and Plemmons (1994, p. 135,137), we directly 
establish NP-hardness of the problems of testing if there is a convex com- 
bination of a given set of matrices, whose real eigenvalues are all positive 
(Dig), or whose inverse is a nonnegative matrix (N38). 

Below we briefly describe three problems whose NP-hardness follows 
from theorem 1, and which are of special interest for applications. 

Problem 2 (Polytopic instability). Given a set of matrices, is there a con- 
vex combination of those that is Hurwitz stable? 

If B,,- is a nonsingular M-matrix, then — B,r is Hurwitz. Therefore, our 
result establishes NP-hardness of the problem of testing the instability of a 
linear time-invariant system with polytopic uncertainty, providing evidence 
for the difficulty of obtaining low-order necessary conditions for this problem 
(Chesi, 2010). The problem of testing the stability of a polytopic uncertain 
system is also NP-hard (Gurvits and Olshevsky, 2009). We refer to Blondel 
and Tsitsiklis (1997) for other complexity results in control theory. 

Problem 3 (Positive determinant). Given a set of symmetric matrices, is 
there a convex combination of those whose determinant is positive? 

NP-hardness of this problem for non-symmetric matrices is established in 
the course of the proof of theorem 1. NP-hardness of the case of symmetric 
matrices follows from a result by Grenet et al. (2011, theorem 5) that shows 
that, for any given n x n matrix X = [xij], one can construct a symmetric 
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matrix Y of dimension 0(n 3 ) with entries in {xjj}U{0; 1; —1; 1/2} such that 
det(X) = det(Y). Note that the problem becomes convex when we constrain 
the feasible region to positive definite matrices (Boyd and Vandenberghe, 
2004, p. 73). 

Problem 4 (Minimize spectral radius). Given a set of nonnegative matri- 
ces, is there a convex combination of those whose spectral radius is less than 
one? 

NP-hardness of this problem follows from theorem 1. Note that the 
matrix B,r in (2) can be written as B,,- = I — M^, where the matrix M,,- = 
Ya=i ^i^i i s a convex combination of nonnegative matrices 



N 



e; + a 
el 1 - i 



(5) 



From Berman and Plemmons (1994, p. 133, definition 1.2) it then follows 
that B,,- is a nonsingular M-matrix if and only if the spectral radius of M,,- 
is strictly less than one. 

This result establishes that the problem of globally minimizing spectral 
radius over a matrix polytope is NP-hard. The spectral radius of a non- 
symmetric real matrix affine function such as NL n is generally a nonconvex 
function of the parameter vector 7r, and the problem of globally optimizing 
the spectral radius is a difficult one (Overton and Womersley, 1988; Han 
et al., 1999; Neumann and Sze, 2007). Our result provides evidence for this 
'difficulty'. NP-hardness of spectral radius minimization was also recently 
established for a class of nonnegative matrix affine functions (Fercoq, 2011). 

We note that the problem becomes convex, and often admits tractable 
exact global optimization algorithms, in the case of symmetric matrices 
(Overton and Womersley, 1988; Fiedler, 1992), or in the case of irreducible 
nonnegative matrices whose entries are posynomial functions of a parameter 
vector 7r. In the latter case, minimizing spectral radius under posynomial 
constraints in 7r can be solved in polynomial time via the Collatz-Wielandt 
formula (Friedland, 1981) and geometric programming (Boyd and Vanden- 
berghe, 2004, p. 165). Note that in our polytope- M-matrix problem the 
simplex equality constraint Y^i=i = 1 is not a posynomial constraint (it 
would be if the equality were replaced by the inequality '<'). 
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